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Abstract: To predict fractal decoded image quality more efficiently, an effective decoded image quality 
prediction method was proposed in this study. In fractal encoding process, the dynamic range of the linear 
correlation coefficients (LCCs) between range blocks and their best-matched domain blocks was greatly 
extended by several outliers which increased uncertainty and resulted in reduced prediction accuracy. To 
remove the interference of outliers, we introduced the effective minimum and maximum of LCCs, which 
provided the effective bottom and top limits of the actual percentage of accumulated collage error 
(EBL-APACE and ETL-APACE), respectively. Further, when EBL-APACE reached a large percentage, the 
average collage error (ACER) can be estimated, and the decoded image quality can be predicted directly. 
Experimental results show that compared with the previous method, the proposed method can provide 
higher prediction accuracy with fewer computations. 
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1. Introduction 


With the rapid development of modern society, people are gradually exposed to more and 
more images in their daily life, which brings great pressure to the transmission and storage of 
images. Efficient image compression technology is one of the important ways to solve the 
above problems. Different from conventional image compression techniques, fractal image 
coding has the characteristics of novel idea, potential high compression ratio, fast decoding, 
and resolution independence, which have attracted the attention of researchers worldwidel!“. 
However, fractal image coding has the drawback of high computational complexity in the 


encoding process, which seriously slows down its practical application. Thus, fast fractal 


image encoding becomes one of the key research directions. In recent years, many fast fractal 
encoding methods were proposed which can be mainly divided into two categories: 1) Local 
block-matching based fast fractal image encoding!*~’!. This type of methods accelerated the 
encoding process by excluding a number of domain blocks and converting global matching 
into local matching. 2) No-search fractal image coding!*!"|. This type of methods assigned the 
best-matched domain blocks directly without searching operations, and real-time encoding 
and higher compression ratio can be realized at the expense of decoded image quality. In 
addition to image compression, fractal image coding has been gradually applied to other 


image processing applications, such as image denoising!!?!®! , image hashing!!?°l, image 
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retrievall?!-31], watermarking’*?°), image super resolution!?”*"!, head pose estimation! 
and MRI image processing!4). 

Predictable quality of decoded images is a novel property of fractal image coding. In 
previous work, the quality of decoded images was predicted only with partial encoding 
process and without decoding process!**!. In this study, an effective method was proposed to 
predict the quality of decoded images more accurately and quickly. First, we found that the 
dynamic range of the linear correlation coefficients between range blocks and their 
best-matched domain blocks was extended greatly by several outliers. Because the number of 
outliers was very small, they contributed little to the accumulated collage error (ACE) of all 
range blocks. Thus, we proposed a novel method to calculate the effective minimum and 
maximum of linear correlation coefficients to remove the interference of outliers, which 
provided the effective bottom and top limits of the actual percentage of accumulated collage 
error (EBL-APACE and ETL-APACE), respectively. Moreover, when EBL-APACE reached 
a large percentage (Such as 90%), EBL-APACE and ETL-APACE determined a small 
dynamic range for the actual percentage of accumulated collage error (APACE), and the 
average collage error (ACER) can be estimated approximately. Further, based on the 
logarithmic relationship between ACER and the fractal decoded image quality, the quality of 
decoded images can be predicted directly. Finally, we adopted four fractal encoding methods 
to assess the performance of the proposed method. Experimental results show that compared 
with the previous method, the proposed method can provide higher prediction accuracy with 


fewer computations. 


This paper is organized as follows: Conventional fractal image coding is reviewed in 
Section 2. The principle of the proposed method is described in Section 3. In Section 4, the 
experimental results are presented and analyzed in detail. The final conclusion is given in 


Section 5. 


2. Conventional fractal image coding 


Fractal image coding aims to establish an iterated function system whose fixed point can 
approach the original image approximately. The first block-based fractal encoding method 
was proposed by Jacquin!'!. First, the original M x N image f was uniformly divided into 
nonoverlapping B x B range blocks Ri, i=1,2,3,....NumR, where NumR denotes the total 
number of range blocks. Then, domain blocks Dj, j=1,2,3,....NumD, can be obtained by 
sliding a 2B x 2B window over the original image from left to right and from top to bottom, 
where NumD denotes the total number of domain blocks. Further, domain blocks were 
contracted to the same size of range blocks, and an extended domain block pool (EDBP) was 
obtained by performing eight isometric transformations on domain blocks. For a given image, 
Figure 1 illustrates the eight isometric transformations which contain rotation through 0 
degree, 90 degree, 180 degree, 270 degree, reflection about the middle horizontal axis, 
reflection about the middle vertical axis, reflection about the main diagonal and reflection 
about the secondary diagonal in order. Finally, the best-matched domain block can be 


obtained by minimizing the following function 


7 ool 
CE(R,) = min min 2 


R,-s,D, -onf ) 
(1) 

i=1,2,3,..., NumR, j =1,2,3,..., NumDx8 
where CE(R;) denotes the collage error of R;. I denotes a B x B matrix whose elements are all 


ones. s; and o; denote scaling and offset coefficients, respectively. D; is the best-matched 


domain block of R;. With the least square method, s; and o; can be calculated as follows 


ead (2) 


i Bran] 


s,=(R,-7, D,-d,1)/|D,-d,1| 


where (*.*) denotes the inner product. 7 and d, are the averages of R;and Dj, respectively. 
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Fig.1 Eight isometric transformations. 


In decoding process, arbitrary M x N image can be selected as the initial image. First, 
the same transformations as those in encoding process were conducted from the best-matched 
domain blocks to corresponding range blocks. One iteration can be completed after all range 
blocks were reconstructed. After about 10 iterations, the decoded image can be obtained. 
Figure 2 (a) selects a rotated 256 x 256 Bridge as the initial image. When we encoded Bridge 
in the encoding process, Figure 3 (b) to (f) illustrate the first five iteration images in decoding 


process. We see that the iteration images converge to Bridge gradually. 


re er À ; 


(a) Initial Image (b) First iteration (c) Second iteration 
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(d) Third iteration (e) Forth iteration (f) Fifth iteration 


Fig.2 Fractal decoding process. 


3. Proposed method 


3.1 Effective minimum and maximum of the linear correlation coefficients between range 
blocks and their best-matched domain blocks 


For Rj, substituting (2) back to (1) yields 


1 2 
CE(R,) = alk -s,D, -04| 


-h(l -at 512, 
wae |x, -7a |p, -a1 
=—|R -71f (I-Loc?) 


where LCC; denotes the linear correlation coefficient between R; and Dj. Because LCC? lies 


between 0 and 1, Eq. (3) satisfies 


j 


R,-FI 


CE(R,) < = (4) 
In Eq. (4), the variance of R; provides the top limit of CE(R;), and thus only the range blocks 
with large variances may result in large collage errors. By Eq. (3), ACE of all range blocks 
can be calculated as 
1 Num 
2 


ACE“ Sen » |R -71 
i=l B i=l 


P[1-(Lcc,y"| (5) 


If we use LCCmax and LCCmin to denote the maximum and minimum of LCC,, 


i=1,2,3,....NumR, respectively, and by Eq. (5), we have 


ACES, < ACE*" < ACES) (6) 
where 
ACEM = Y CER) = Y |R -7| [1-(LCCu) | 
Min L i / Min B? = i i Max 
= 2) 
ACE vin = 2, CE(R; ) ax z B 2, |R, pa all [1- (LCC, } | 


Linear correlation coefficients 
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Fig.3 Distribution of linear correlation coefficients. 


In encoding process, we sorted the range blocks by their variances from largest to smallest 


and encoded them in order. For Bridge, Figure 3 shows the distribution of LCCs. In Eq. (7), 


All 


LCCmax and LCCmin provide the bottom and top limits of ACE“!, ACE,;, and ACE\,.., 


J respectively. However, we observe that the dynamic range of LCCs is mainly determined by 
several outliers which greatly extend the dynamic range of LCCs. In Fig.3, all the LCCs 


below the red line, i.e., smaller than 0.6, are labeled by “o”. The outliers make LCCmin extend 


to below 0.6. Because the outliers’ ACE only comprise 1.47% of total collage errors, their 
effect can be ignored. To remove the interference of outliers and effectively estimate the 
dynamic range of ACE“", we adopt the following method to compute the effective minimum 
and maximum of LCCs: 

For all range blocks, we established a histogram of collage errors with respect to LCCs. 
If we had NumC range blocks Ry, k=1,2,3,....NumC, whose LCCs were all Iccs. We used 


P(lcc) to represent the normalized version of the sum of the collage errors with the same LCC 
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as 


NumC 
Bile = > CE, (ce) /ACE™ Icc € [0,1] 


k=1 


where CEx(Icc) denotes the collage error of Rx with LCC=Icc. From Eq. (8), we know that 


P(Icc) actually represents the probability density function o 


respect to LCCs. 


f collage errors (PDF-CE) with 


Probability density function of collage errors 
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Fig.4 PDF-CE versus LCCs. 
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(a) Fitting result of the left part and its reflection.(b) Fitting result of the left part and its reflection. 


Fig.5 Fitting results for the left and right parts of PDF-CE in Fig.5, respectively. 
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Fig.6 Fitting results of PDF-CE. 


Figure 4 shows PDF-CE for Bridge which is asymmetric. We plotted a vertical line (VL) 
at lec=7’ where PDE-CE reached its peak value. VL divides PDF-CE into the left and right 


parts which comprise Tret% and Trig% of ACE4", respectively, and we have 


Tien YtT, 


Right %= 1 00% (9) 
We took out the left part and reflect it about VL, a symmetric curve in Fig.5 (a) was 
obtained. Further, by reflecting the right part about VL, we also obtained a symmetric curve 


in Fig.5 (b). Then, we adopted the following two Gaussian functions, GFret(lcc) and 


GFRignt(lcc), to fit the two curves in Fig.5 (a) and (b), respectively. 


(10) 


where a3 and f3 represent the standard deviations of the curves in Fig.5 (a) and (b), 
respectively. a2 and f2 are both equal to T. We know that Gaussian distributions concentrate 
95.5% x 2TLet% and 95.5% X 2Trign% of ACE^! in the intervals of [a2-2a3, a2+2a3] and 
[Po-283, B2t2B3] in Fig.5 (a) and (b), respectively. Thus, [a2-2a3, a2] and [fo, B2+2f3] 
concentrate 95.5% x Tren% and 95.5% Xx Trig% of ACE", respectively. Further, we took out 
the left and right parts of the fitting curve in Fig.5 (a) and (b), respectively, and combined 
them together in Fig.6. By Eq. (9), we know that the interval [a2-2a3, /2+2/3] concentrate 
95.5% of ACE". Finally, we define a2-2a3 and f2+2f3 as the effective minimum and 


maximum of LCCs, LCCmin Effective and LCCMax Effective, respectively, as 


eae = a, -20 


11 
LCC Max_Effective = LA i 2 B; ( ) 


LCCmin Effective and LCCmax Effective can effectively remove the effect of outliers. Because the 
the dynamic range of LCCs is narrowed, we also have the inequalities as 


| LCC Min Effective > LCC yin ( 1 2) 


LCC < LECC 


Max Effective Max 


3.2 Effective bottom and top limits of accumulated collage error 


With LCCmin Effective and LCCmax Effective, Eq. (7) can be modified as 


NumR 1 


NumR 
ACE in = z CEC(R;) min Effective = B > |R, = Fal E = e Ji 
i=l i=l (13) 


NumR 1 


NumR 
All = = 
ACE a Effective ~~ Ds CE(R,) max. Effective a B? T | 
i=l i=l 


R-7I 


2 2 
| f E e ) | 
We compared ACE! with ACE®in pnecive as 


1 NumR 
ACE =~ ACE ae = B? 2 | CE(R,) G CE(R; ) Min_Betectve | 
i=l 
(14) 


1 NumR 


ay) 


i=l 


R, Ez Fl e i pP (LCC, Y | 


The range blocks can be divided into two categories: Range blocks whose LCCs are above 
and below LCCMax Effective, respectively. We use Num*ve and Num®e!’’ to denote the numbers 
of range blocks whose LCCs are above and below LCCmax Effective, respectively. Then, Eq. (14) 


can be rewritten as 


NumAbove 
ACE” a AGE rere a = |R, Fy rI | eee oe j z; (LCC, y | a 
i=l 
aes oa (15) 
2 
=o > | R, ae rI i pae z (LCC, | 
i=l 

>0 


Only the LCCs of several outliers satisfy LCC>LCCwax Effective, and the others satisfy 


LCC;<LCCmMax Effective- For Bridge, Num4?ove=39, Num®*w=4.057, and 
ACE! — ACE Win tttective = 529,480 > 0. Generally, the number of LCCs below LCCMax Extective 


is more than ten times larger than above LCCmax Effective. Thus, in Eq. (15), in the case of 


All All 
Num®ow>>Num4?ove, ACE“ — ACE isin effective > 9, and we have 


ACE“ a AGH pages ( 1 6) 


Then, we compared ACE^! with ACE Six grtective aS 


1 NumR 


ACE ie Gages ~ACE™=—5 > [CECR e y CE(R,) | 
i=l 
1 NumR = 
a 2, |R, = rd 


s ; À (17) 
ee ) E (ECC 2s mee ) | 
Similar with Eq. (15), we divided range blocks into two categories: The range blocks whose 


LCCs are above and below LCCnin Fitting, respectively. We use Num*%ve and Num®ow to 
denote the numbers of range blocks whose LCCs are above and below LCCmin Effective, 


respectively. Then, Eq. (17) can be rewritten as 


1 Num^’o”e 


ACE aee ry ACE*!= B? 2 |R, B rI keo y ~ (DEC aero y | 


i=l 


>0 
1 Num Bebow (1 8) 


B? 2, |R, -nI ; iee I) F (EC Y | 


<0 


Only the LCCs of several outliers satisfy LCC<LCCnin Effective, and the others satisfy 


LCC>LCCmin Effective, For Bridge, Num/eve=4,026, NumPeow=70, and 


ACH emai —ACE*“" = 760,276 > 0. Thus, in Eq. (18), in the case of Num*%ve>>NumBPeow, 


we have 
ACB Geen z ACE®™ (19) 
By Eqs. (16) and (19), we have 
2 
ACE a = ACE*" < ACE o ( 1 


Thus, the effective minimum and maximum of LCCs, LCCmin Effective and LCCMax Effective, can 
provide the effective top and bottom limits of ACE“", ACE m and ACE ena. 3 
respectively. 


3.3 Decoded image quality prediction 


In encoding process, range blocks were divided into two categories: Coded range blocks and 


uncoded range blocks. For coded ones, if the number of them is large enough, we can use the 


effective minimum and maximum of the LCCs of coded range blocks, LCC" t roire and 

LCC yin pttecive > t0 approximate those of all range blocks, LCCmin Etrective and LCCmax Ftfective, aS 
Coded Coded 

LOC hin Effective oe LCC Min Effective > LOC iiex Effective fee LOC wren Eitective (2 1) 


Moreover, the effective minimum and maximum of the LCCs of uncoded range blocks, 


Uncoded Uncoded 
LCC hin recive aNd LCCh Effective , can be represented by those of all range blocks as 


10 


Lees -LCC 


Min Effective ? Max_Effective 


ECEN =LCC 


Min Effective 


(22) 


Max_Effective 
By Eqs. (21) and (22), we can use EOC and CC aes to approximate 


LCC Uncoded and LC C Uncoded as 


Min Effective Max_Effective 


ECCS x LCC Urceded Leese ~ LECE (23) 


Min Effective Min Effective ? Max_Effective Max Effective 


Further, by Eq. (3), the collage error of arbitrary uncoded range block can be denoted as 


CE(R™) = jeee -71| E roe | i=1,2,3,..., Non (24) 


where Num”*°°e¢ denotes the number of uncoded range blocks. LCC} denotes the LCC 


of R™** Then, the ACE of uncoded range blocks can be described as 


Num Uncoded 


ACE Uncoded = > CE(R Vas) = 


i=l 


1 NumUncoded ; y 
aa pa Roe -71| parece) | (25) 
i=l 


We assume that the result in Subsection 3.2 is also suitable for the uncoded range blocks, 1.e., 


the effective minimum and maximum of the LCCs of uncoded range blocks, [Rel eras 
and LCCyin'Snective » Can provide the effective top and bottom limits for ACEY@4 as 


ACE Uncoded < ACE Uncoded < ACE Uncoded (26) 


Min Effective Max_Effective 


where 


Uncoded Uncoded 


Num 1 Num 2 2 
Uncoded as Uncoded a. Uncoded Uncoded 
ACE Min Effective = > CE( R; ) Min Effective z B? p2 |k; z rI | i p Ee ) | 
i=l 
Uncoded 


i=l 


Num Num Urcoded 


Uncoded A Uncoded LAD 1 Un coded — yll? Uncoded 2 
ACE Max Effective a > CE( R, ) Max ea g B? 2 |z; = 71| ii J E ey. ) 


i=l i=l 


(27) 
By Eq. (23), Eq. (27) can be rewritten as 


Uncoded —_ oe Uncoded oF 1 a Uncoded —,||2 Coded 2 
ACE Min Effective ~ 2, CE(R; ) Min Effective = B 2, |R, =< r I | i z (cere ) | 

Uncoded yee Uncoded 1 i Uncoded = 2 Coded 2 
ACE wax Effective 5 2, CE(R; ) Max Effective z B 2, |R; T 71| i F (CC a aa ) | 


(28) 
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Moreover, if the actual percentage of accumulated collage error (APACE) is defined as the 


ratio of ACE©*4 to ACE“! as 


ACE Coded ACE Coded 
APACE= = 29 
ACE*" ACE&“#4 m ACE Uncoded ( ) 


By Eqs. (26) and (29), we obtain the effective bottom and top limits of the actual percentage 
of accumulated collage error, EBL-APACE and ETL-APACE, as 


ACES : ACEH 2 ACE Coded 
ACE Sa an ACE Nee Ş ACE Bn an ACE aay ACE pp a ACE n T (3 0) 
EBL-APACE APACE ETL-APACE 


From Eq. (30), we know that ETL-APACE<1 and EBL-APACE<ETL-APACE. When 
EBL-APACE reaches a large percentage (Such as 90%), EBL-APACE can approach 
ETL-APACE sufficiently, which implies that APACE lies in a rather small dynamic range and 


can be approximately estimated as 


EBL-APACE + ETL - APACE 
9 (31) 


Moreover, the average collage error (ACER) of all range blocks can be computed as 


APACE Estimated _ 


ACE a 
ACER = ees >, CE(R,) /NumR (32) 
u 


i=l 


By Eqs. (32), (29), and (31), we can estimate ACER as 


ACE Coded /APACE Estimated GB 3) 
NumR 


ACER 7% 


The numerator of Eq. (33) represents the estimated version of ACE". In addition, the peak 
signal-to-noise ratio (PSNR) is used to measure the decoded image quality as 
ane Decoded 2 
PSNR = log, 4 255° Fae rA CEN ae CD) (34) 
Mx N i= 
where f and fPeoded represent the input and decoded images, respectively. Moreover, there 
exists a logarithmic relationship between ACER and the PSNR of the decoded image as!°°! 


PSNR = £, + f,log,, (ACER) (35) 


where fı and f2 are both constant values. Finally, by Eqs. (33) and (35), we can directly 
predict the PSNR of the decoded image. 


12 
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4. Experiments 


We adopted six 256 x 256 images (in Fig.7), Peppers, Airplane, Boat, Bridge, House, and 
Lake, in our experiment. Four fractal image coding methods were used to assess the 
performance of the proposed method!'»-7). The scaling and offset coefficients, s and o, were 
quantized by 5 and 7 bits, respectively. The procedures of the proposed method are listed as 
follows: 

Step 1: Divide the input image into range and domain blocks, and sort the range blocks 
by their variances from largest to smallest. 

Step 2: Take out one uncoded range block and encode it. If EBL-APACE is smaller than 
90%, turn back to Step 2. If not, calculate ETL-APACE and go to Step 3. 


Step 3: Calculate ACER by Eq. (33), and predict the decoded image quality by Eq. (35). 
R E3 | eed a | | 
m 


(a)Peppers (b) Airplane (c) Boat 


(d) Bridge (e) House (f) Lake 
Fig.7 Six test images. 
When the size of range blocks is set to be 4 x 4 and 8 x 8, respectively, Table 1, 2, 3, and 
4 list the experimental results of Jacquin’s, Chaurasia’s, Zheng’s, and Gupta’s methods, 


respectively!'*-7!, In Tab.1, when the size of range blocks is 4 x 4, for Jacquin’s method, the 


13 


second row shows the PSNRs of decoded images. In the third row, we encode all the range 
blocks, conduct the block-matching 4,096 times for all test images, and the percentage of 
computations is considered as 100%. For the previous method, the fourth to sixth rows list 
the predicted PSNRs, deviations between the predicted and actual PSNRs, and the percentage 
of computations with respect to Jacquin’s method, respectively. Similarly, the seventh to 
ninth rows show the counterparts for the proposed method. From Tab.1, we see that if the size 
of range blocks is 4 x 4, the average deviations of the previous and proposed methods are 
0.08dB and 0.05dB, respectively, the average percentages of computations are 52.99% and 
47.76%, respectively. If the size of range blocks is 8 x 8, the average deviations of the 
previous and proposed methods are 0.10dB and 0.09dB, respectively, the percentage of 
computations needed are 52.62% and 50.73%, respectively. We see that compared with the 


previous method, the proposed method can provide higher prediction accuracy with fewer 


computations. The analyses are as follows: In previous method, LCCH and LCCHES 


were used to approximate LCCH and LCCH! and provided the bottom and top 
limits of APACE as 
AGE == NCE 
APACE = APACE, = 
Min ACEC + ACE Yreoded > Max ACEC! + ACE Yed (36) 
where 
Uncoded __ e Uncoded a l i Uncoded = 2 Coded 2] 
ACEN = CER a QL RP FA] 1- (LCC) 
N aia 1 N a (37) 
ncode i ncode i ncode e oded \2 ] 
ACE Di >a CE(R? 3 ee z B |R -c 71| i C) 
i=l i=l 4 


Similar with Eq. (12), for coded range blocks, we also have 


LGC ens > TEN (3 8) 
LCG o es < ECC 
Then, by comparing Eq. (28) and (37), we have 
APACE, ,,,<EBL-APACE , ETL-APACE<APACE,,,.. (39) 


By Eq. (39), APACEmin is always smaller than EBL-APACE. When EBL-APACE reaches 


14 


90% in the proposed method, APACEmin is smaller than 90%, which implies that in the 
previous method, we need to encode extra range blocks to make APACEmin reach 90%. Thus, 
a number of computations in the proposed method are saved. When the size of range blocks 
is set at 4 x 4, we totally have 4,096 range blocks. For Peppers, in the previous method, when 
APACE min reaches 90%, we need to encode 2,327 range blocks which are represented by the 
red and white boxes in Fig.8 (a) and comprise 56.81% of total range blocks. In the proposed 
method, when EBL-APACE reaches 90%, we only need to encode 2,002 range blocks which 
are represented by white boxes in Fig.8 (a) and comprise 48.88% of total range blocks. This 
implies that the computations of 325 range blocks which are represented by red boxes are 
saved. When the size of range blocks is set at 8 x 8, we totally have 1,024 range blocks. In 
the previous method, when APACEmin reaches 90%, we need to encode 529 range blocks 
which are represented by the red and white boxes in Fig.9 (a) and comprise 51.66% of total 
range blocks. In the proposed method, when EBL-APACE reaches 90%, we need to encode 
487 range blocks which are represented by the white boxes in Fig.9 (a) and comprise 47.56% 
of total range blocks. This implies that the computations of 42 range blocks which are 
represented by red boxes are saved. Moreover, By Eq. (39), we also know that compared with 
APACEmin to APACEmax, EBL-APACE to ETL-APACE provides smaller dynamic range for 
estimating APACE, and the proposed method is expected to provide higher prediction 
accuracy. 

In Tab.2, for Chaurasia’s method, if the size of range blocks is 4 x 4, the average 
deviations of the previous and proposed methods are both 0.07dB and 0.05dB, respectively. 
The average percentages of computations are 54.14% and 48.12%, respectively. If the size of 
range blocks is 8 x 8, the average deviations of the previous and proposed methods are both 
0.17dB. The average percentages of computations are 52.44% and 50.31%, respectively. In 
Tab.3, for Zheng’s method, if the size of range blocks is 4 x 4, the average deviations of the 
previous and proposed methods are 0.14dB and 0.10dB, respectively. The average 
percentages of computations are 53.87% and 47.70%, respectively. If the size of range blocks 
is 8 x 8, the average deviations of the previous and proposed methods are 0.11dB and 0.09dB, 
respectively, The average percentages of computations are 52.60% and 50.68%, respectively. 
In Tab.4, for Gupta’s method, if the size of range blocks is 4 x 4, the average deviations of the 
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previous and proposed methods are 0.08dB and 0.04dB, respectively. The average 
percentages of computations are 53.24% and 48.13%, respectively. If the size of range blocks 
is 8 x 8, the average deviations of the previous and proposed methods are 0.15dB and 0.14dB, 
respectively. The average percentages of computations are 52.80% and 50.90%, respectively. 
In summary, compared with the previous method, the proposed method can predict the 


decoded image quality more accurately with fewer computations. 


(a)Peppers (b) Airplane (c) Boat 


(d)Bridge (e) House (f) Lake 
Fig.8 When the size of range blocks is 4 x 4, white and red boxes represent the range blocks encoded in the 


previous method, and white boxes represent the range blocks encoded in the proposed method. 
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(a)Peppers (b) Airplane (c) Boat 


16 


202308.00720v1 


ChinaXiv 


(d)Bridge 


(e) House 


Fig.9 When the size of range blocks is 8 x 8, white and red boxes represent the range blocks encoded in the 


previous method, and white boxes represent the range blocks encoded in the proposed method. 


Tab.1 Performance comparison between the previous and proposed methods for Jacquin’s method", 


Test images Peppers Airplane Boat Bridge House Lake Average 
R Decoded(dB) 31.80 29.46 28.28 25.44 27.56 D7: x 
Jacquin’s!"] i 
Computations(%) 100 100 100 100 100 100 x 
Predicted(dB) 31.64 29.39 28.30 2553 27.52 27.64 x 
Previous! Deviation(dB) 0.16 0.07 0.02 0.11 0.04 0.08 0.08 
Computations(%) 56.81 35.38 53.93 71.17 50.27 50.37 52.99 
Predicted(dB) 31.71 29.42 28.32 25.56 27.55 VIKIA x 
Proposed Deviation(dB) 0.09 0.04 0.04 0.12 0.01 0.01 0.05 
Computations(%) 48.88 32.50 48.24 67.11 46.17 43.65 47.76 
' Decoded(dB) 25.28 23.17 23.28 21.19 22.37 22.18 x 
Jacquin’s!") ; 
Computations(%) 100 100 100 100 100 100 x 
Predicted(dB) 25.27 23.14 23.17 21.11 22°52 21.93 x 
Previous)! Deviation(dB) 0.01 0.03 0.11 0.08 0.15 0.25 0.11 
Computations(%) 51.66 36.91 53.71 72.17 51.37 49.90 52.62 
Predicted(dB) 25.28 23.11 23.16 21.07 22.48 2192 x 
Proposed Deviation(dB) 0.00 0.06 0.12 0.12 0.11 0.16 0.10 
Computations(%) 47.56 35.94 51.95 70.80 50.00 48.14 50.73 


Tab.2 Performance comparison between the previous and proposed methods for Chaurasia’s method", 


Test images Peppers Airplane Boat Bridge House Lake Average 

Chaurasia’s Decoded(dB) 31.47 29.37 28.04 25.16 27.08 27.28 x 
[5] Computations(%) 100 100 100 100 100 100 x 
Predicted(dB) 31.30 29.23 28.00 25.12 27.07 27.26 x 

Previous)! Deviation(dB) 0.17 0.14 0.04 0.04 0.01 0.02 0.07 

Computations(%) 59.72 35.25 55.98 72.02 50.10 51.76 54.14 
Predicted(dB) 31.40 29.26 28.05 25.14 27.09 27.34 x 

Proposed Deviation(dB) 0.07 0.11 0.01 0.02 0.01 0.06 0.05 

Computations(%) 49.37 32.74 48.75 67.77 46.17 43.92 48.12 
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Chaurasia’s Decoded(dB) 24.97 23.15 23.08 20.64 22.19 21.82 x 
[5] Computations(%) 100 100 100 100 100 100 x 
Predicted(dB) 24.85 23.03 22.89 20.87 2237 21.64 x 
Previous’ Deviation(dB) 0.12 0.12 0.19 0.23 0.18 0.18 0.17 
Computations(%) 50.39 36.91 54.30 72.07 51.07 49.90 52.44 
Predicted(dB) 24.86 23.00 22.90 20.85 22.33 21.61 x 
Proposed Deviation(dB) 0.11 0.15 0.18 0.21 0.14 0.21 0.17 
Computations(%) 45.90 35.74 51.66 70.61 50.00 47.95 50.31 
Tab.3 Performance comparison between the previous and proposed methods for Zheng’s method", 
Test images Peppers Airplane Boat Bridge House Lake Average 
Decoded(dB) 31.78 29.29 28.04 25.06 27.46 27.70 x 
Zheng’s!* : 
Computations(%) 100 100 100 100 100 100 x 
Predicted(dB) 31.61 29.19 28.07 25.43 27.43 27.56 x 
Previous*] Deviation(dB) 0.17 0.10 0.03 0.37 0.03 0.14 0.14 
Computations(%) 57.50 35.99 55.08 70.78 50.02 53.86 53.87 
Predicted(dB) 31.69 29.24 28.10 25.44 27.46 27.67 x 
Proposed Deviation(dB) 0.09 0.05 0.06 0.38 0.00 0.03 0.10 
Computations(%) 48.83 32.37 48.90 66.65 45.92 43.55 47.70 
Decoded(dB) 25.24 23.06 23:25 20.99 22.36 PARTE x 
Zheng’s'*] i 
Computations(%) 100 100 100 100 100 100 x 
Predicted(dB) 25.26 23.14 23.14 21.11 225) 21.93 x 
Previous!?>] Deviation(dB) 0.02 0.08 0.11 0.12 0.15 0.16 0.11 
Computations(%) 51.66 36.91 53.61 72.17 51.37 49.90 52.60 
Predicted(dB) 25.28 23.11 23.13 21.07 22.47 21791 x 
Proposed Deviation(dB) 0.04 0.05 0.12 0.08 0.11 0.14 0.09 
Computations(%) 47.56 35.84 51.76 70.80 50.00 48.14 50.68 
Tab.4 Performance comparison between the previous and proposed methods for Gupta’s method”. 
Test images Peppers Airplane Boat Bridge House Lake Average 
Decoded(dB) 31.75 29.29 28.16 25.32 27.48 27.66 x 
Gupta’s!7] i 
Computations(%) 100 100 100 100 100 100 x 
Predicted(dB) 31.57 29.26 28.14 25.38 27.41 27.56 x 
Previous“ Deviation(dB) 0.18 0.03 0.02 0.06 0.07 0.10 0.08 
Computations(%) 58.20 35.69 54.39 71.41 50.66 49.07 53.24 
Predicted(dB) 31.66 29.29 28.15 25.38 27.43 27.61 x 
Proposed Deviation(dB) 0.09 0.00 0.01 0.06 0.05 0.05 0.04 
Computations(%) 48.58 32.81 49.71 67.29 46.53 43.85 48.13 
Decoded(dB) 25.26 23.08 2322 20.88 22.20 22.13 x 
Gupta’s!7] j 
Computations(%) 100 100 100 100 100 100 x 
Predicted(dB) 25.19 23.00 23.06 21.00 22.39 21.85 x 
Previous>] Deviation(dB) 0.07 0.08 0.16 0.12 0.19 0.28 0.15 
Computations(%) 51.56 37.30 54.20 72.36 51.17 50.20 52.80 
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Predicted(dB) 25.19 22.97 23.06 20.96 22.34 21.84 x 
Proposed Deviation(dB) 0.07 0.11 0.16 0.08 0.14 0.29 0.14 
Computations(%) 48.34 36.04 51.76 70.70 50.10 48.44 50.90 


5. Conclusion 


An effective method for predicting the fractal decoded image quality was proposed in this 
paper. First, we introduced the effective minimum and maximum of LCCs, which provided 
EBL-APACE and ETL-APACE for APACE. Then, decoded image quality can be directly 
predicted while EBL-APACE reached a large percentage. Experimental results show that 
compared with the previous method, the proposed method can provide higher prediction 
accuracy with fewer computations. In future work, we will focus our attention on the 
following two aspects: 

a) Maintain prediction accuracy. By removing the interference of the outliers of LCCs, we 
have obtained satisfying prediction accuracy. In future work, we will try to maintain the 
current prediction accuracy. 

b) Shorten prediction process. Compared with the previous method, the number of 
computations in the proposed method has been reduced. In future work, we will try to reduce 


the number of computations further. 
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